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Abstract. Recent progress in the theory of solar and stellar dynamos is reviewed. Particular emphasis is placed on the 
mean-field theory which tries to describe the collective behavior of the magnetic field. In order to understand solar and 
stellar activity, a quantitatively reliable theory is necessary. Much of the new developments center around magnetic helicity 
conservation which is seen to be important in numerical simulations. Only a dynamical, explicitly time dependent theory of 
a-quenching is able to describe this behavior correctly. 
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1. Introduction 

Starspot activity is presumably driven by some kind of dy- 
namo process. Many stars show magnetic field patterns ex- 
tending over scales of up to 30° in diameter. The commonly 
used tool to model such magnetic activity is the mean-field 
dynamo. Although mean-field theory has been used over sev- 
eral decades there have recently been substantial develop- 
ments concerning the basic nonlinearity of dynamo theory. 
It is the purpose of this review to highlight these recent de- 
velopments in the light of applications to stars. 

2. Stellar dynamos: spots and cycles 

We usually think of star spots as rather extended dark and 
strongly magnetized areas on a stellar surface. Observable 
spots are much bigger than sunspots. They may in fact be 
so big that the spots themselves have sometimes been iden- 
tified with solutions of the mean-field dynamo equations 
(e.g., Moss et al. 1995). This is in contrast to the much 
smaller sunspots which instead show collective behavior in 
that sunspot pairs have a systematic orientation and preferen- 
tial location which changes with the solar cycle. 

The working hypothesis is that extended star spots are just 
the extremes of a broad range of possibilities from small to 
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large spots. Stellar parameters can change over a consider- 
able range and there is scope that different types of behavior 
can be identified with different solutions of the mean-field 
dynamo equations. Very exiting is the possibility of nonax- 
isymmetric solutions, possibly with cyclic nonmigratory al- 
ternations of their polarity (the so-called flip-flop effect, see 
Jetsu et al. 1999 and references therein). 

Already among the more solar-like stars there is a lot to 
be learned about the dependence of the period of the activ- 
ity cycle on rotation rate and spectral type (cf. Baliunas et al. 
1995). An interesting possibility is the suggestion that stellar 
activity behavior may change with age (Brandenburg, Saar & 
Turpin 1998). The very young and more active stars show ex- 
tremely long cycles (3^4 orders of magnitude longer than the 
rotation period) whilst older inactive stars like the sun show 
shorter cycle periods that are just a few hundred times longer 
than the rotation period. These different types of behavior can 
be classified according to their location in the Rossby number 
versus frequency ratio diagram (Saar & Brandenburg 1999). 

In order to make progress in understanding these vari- 
ous possible behaviors it is crucial to work with a reliable 
theory that has predictive power. Mean-field dynamo theory 
has frequently been used as a rather arbitrary theory. Being 
based on some ad hoc assumptions, much of its predictive 
power is questionable. Particular controversy was caused by 
the ill-known contributions of small scale fields which may 
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catastrophically quench the a-effect (Vainshtein & Cattaneo 
1992, Kulsrud & Andersen 1992), which is thought to be re- 
sponsible for driving the large scale field. However, signifi- 
cant advances in recent years are now beginning to shed some 
light on apparently conflicting earlier results on what the fi- 
nal saturation field strength will be. It is likely that progress 
will come about in two stages. In the first stage we will have 
to make sure that mean-field theory works correctly in the 
parameter regime that can be tested using simulations. In 
the second stage we have to extrapolate the theory from the 
regime that is tested numerically to the regime that is of astro- 
physical interest. At the moment we are still struggling with 
the first objective. 



3. Mean-field theory: it's all about quenching 

As far as the selection of different modes of symmetry is 
concerned, there has been some partial success in finding 
agreement between simulations and linear dynamo theory. 
We mention here the results of local simulations of accretion 
discs in a shearing box approximation: changing the upper 
and lower boundary conditions from a normal field ("pseudo- 
vacuum") condition to a perfect conductor condition changes 
the behavior from an oscillatory mode with symmetric field 
about the equator to a non-oscillatory mode with antisym- 
metric field about the equator. The same change of behavior 
is also seen in mean-field models using the same cartesian ge- 
ometry. This result, which has been described in more detail 
in earlier papers (e.g. Brandenburg 1998), lends some sup- 
port to the basic idea of using mean-field theory to describe 
the results of simulations of hydromagnetic turbulence under 
the influence of rotation and stratification and in the presence 
of boundary conditions. 

More serious concern comes from the effects of nonlin- 
earity. Broadly speaking, nonlinearity has to do with strong 
magnetic fields, where the magnetic energy density ap- 
proaches the kinetic energy density of the turbulence. At 
the bottom of the solar convection zone, the correspond- 
ing equipartition field strength is B cq = 4 ... 8 kG. On the ^_ 
one hand, magnetic fields of this strength may actually be 
required for the dynamo to operate. Babcock (1961) and 
Leighton ( 1 969) proposed that magnetic fields of this strength 
will become buoyant and produce, under the influence of the 
Coriolis force, a systematic tilt as flux tubes emerge at the sur- 
face to form a sunspot pair. In many ways magnetic buoyancy 
is similar to thermal buoyancy and both lead to an a-effect 
(Parker 1955, Steenbeck, Krause & Radler 1966). However, 
Piddington (1972) has argued that, when the magnetic field 
approaches the equipartition value, it would be impossible to 
entangle and diffuse the magnetic field. This led later to the 
idea of catastrophic suppression of turbulent magnetic diffu- 
sivity (Cattaneo & Vainshtein 1991) and, by analogy, to the 
proposal of catastrophic suppression of the a-effect (Vain- 
shtein & Cattaneo 1992). Simulations of Tao et al. (1993) 
and Cattaneo & Hughes (1996) show that in the presence of 



an imposed magnetic field, Bo, the a-effect depends on Bo 
like 

= m 
a l + R m Bl/B^ (L) 

where R m is the magnetic Reynolds number which is large 
(10 8 '" 9 for the sun). Thus, for equipartition field strengths, 
Bq ~ B cq , a would be 8 to 9 orders of magnitude below the 
kinematic (unquenched) value ao, i.e. 

a ~ ao-R^ 1 — > as i? m — > oo. (2) 

Over the past ten years there has been an increasing 
amount of activity in trying to understand the value of a in 
the nonlinear regime. Work by Gruzinov & Diamond (1994, 
1995, 1996) and Bhattacharjee & Yuan (1995) has basi- 
cally confirmed Eq. (|l]). Gruzinov & Diamond (1994) did 
find however that the turbulent magnetic diffusivity is only 
quenched in two-dimensional configurations, which was ex- 
actly the case considered numerically by Cattaneo & Vain- 
shtein (1991). Although Gruzinov & Diamond (1994) did 
agree with the conclusion of catastrophic a-quenching, they 
found actually a slightly different form of Eq. (|l|), which can 
be written as 

a + i2 m J7tMoJ • B/fl, 
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eq 



i 



(3) 



where J = V x B//xo is the mean current density and /io 
the vacuum permeability. Obviously, when the mean field is 
spatially uniform, B = Bo = const, then J = and Eqs ([!]) 
and ([3]) agree with each other. 

In real astrophysical bodies a will always be a tensor 
(e.g., Steenbeck et al. 1966, Rudiger & Kitchatinov 1993, Ro- 
gachevskii & Kleeorin 2001). However, much of the work on 
the nonlinear a-effect comes from considering periodic boxes 
where the a tensor can indeed be isotropic (e.g. Field, Black- 
man & Chou 1999). There is a priori no reason to assume 
that the a-effect in a periodic box is different from that in a 
nonperiodic box. Furthermore, periodic boxes are conceptu- 
ally and computationally significantly easier than boxes with 
boundaries. If no large scale field is imposed, helical turbu- 
lence can still drive a large scale field which itself is helical. 
A prototype of such a field is 

B = B (sinfciz, cosfciz, 0), (4) 

where k\ = 2ir/L is the smallest wavenumber in a box of 
size L 3 . Other directions and additional phase shifts are pos- 
sible; see Brandenburg (2001, hereafter B01); see Fig. [j] An- 
imations of x, y, and z slices of the generated magnetic field, 
together with the corresponding power spectra of kinetic and 
magnetic energy, as well as magnetic helicity (normalized by 
k/2) are shown in the attached Movies 1-3. 

A helical field of the form (Q) is called a Beltrami field. 
The current density of such a field no longer vanishes; in fact, 
/UoJ = k%B, so J and B are aligned and 

(UqJ • B = kiH 2 = const. (5) 
Thus, in the large magnetic Reynolds number limit, Eq. (Q) 
becomes 

a = 0{r}%ki) ^ as i? m — > 00. (6) 
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Fig. 1. Visualization of the magnetic field in a three- 
dimensional simulation of helically forced turbulence. The 
turbulent magnetic field is modulated by a slowly varying 
component that is force-free. 



This highlights the great ambiguity in concluding anything 
about a-quenching from oversimplified experiments. [In the 
discussion above we have assumed that a > 0. If a < 0, 
as is the case in B01, then both a and J • B are also negative 
and Eq. (|3|) is unchanged.] 

4. Relation to magnetic helicity 

There is a strong connection between a-quenching and mag- 
netic helicity conservation. Again, we consider the case of a 
periodic box, for which it is easy to show that the magnetic 
helicity 



H = f A BdV = (A-B)y 



(7) 



is perfectly conserved in the limit of infinite magnetic 
Reynolds number. Here, A is the magnetic vector potential, 
so the magnetic field is B = V x A, and angular brackets 
denote volume averages over the full box. If we start with a 
very weak seed magnetic field, the initial magnetic helicity 
must also be small and will therefore always remain small if 
H is conserved. Thus, a large-scale helical field of the form 
(Q) is only compatible with conservation of magnetic helicity 
if there is an equal amount of magnetic helicity of the oppo- 
site sign in the small scales, i.e. 

(A • B) = (A • B) + (a • b) = (early times), (8) 

where A and B have been split up into their mean and fluc- 
tuating components, 



B = B 



(9) 



Overbars refer to the mean-field obtained by horizontal or 
azimuthal averaging, for example. Equation (||) is a crucial 
condition that must be obeyed by any mean-field theory in 



the large R m limit on time scales shorter than the resistive 
time scale. 

To our knowledge, there have been two approaches to in- 
corporate magnetic helicity conservation into mean-field dy- 
namo theory. One is to express the mean turbulent electromo- 
tive force, £ = u x b, as a divergence term (Bhattacharjee & 
Hameiri 1986, see also Boozer 1993) and the other is to mod- 
ify the feedback onto the a-effect such that Eq. (|8|) is satisfied 
on short enough time scales. The latter approach goes back 
to Kleeorin & Ruzmaikin (1982) and Kleeorin et al. (1995), 
and has recently been revived by Field & Blackman (2002), 
Blackman & Brandenburg (2002), and Subramanian (2002). 
In the following we briefly outline the basic idea. 

All we know is that in a closed or periodic domain the 
magnetic helicity evolves according to 

A(A-B)=-2ry M o(J-B), (10) 

where (A • B)V and (J • B)V are magnetic and current he- 
licities, respectively, and V is the volume. At the same time 
we have to have some theory for the evolution of the mean 
magnetic field. The mean-field ail dynamo equations can be 
written in the form 
BB 

— = Vx [UxB+aB-^ + ^oJ], (11) 

From this we can construct the evolution equation for the 
magnetic helicity of the mean field, 

^(A-B)=2(£-B)-2 Wo <J-B), (12) 
where 



£ = aB - ?7tMoJ, 



(13) 



is the mean turbulent electromotive force under the assump- 
tion of isotropy. Subtracting Eq. (|l2|) from Eq. (|io|), we obtain 
the evolution equation for the magnetic helicity of the fluctu- 
ating field, 

A( a .b)=-2(£-B)-2 W o{j-b}. (14) 
at 

This equation has to be solved simultaneously with the usual 
mean field equation. At the moment, however, it is not yet 
fully coupled to the mean field equation. In fact, any kind of 
coupling, for example of the form 

a = a«a-b),(j-b}) (15) 

would suffice. A similar relation could in principle also be 
applied to the turbulent magnetic diffusivity, %. However, in 
contrast to %, (a-b) and (j b) are pseudo-scalars and change 
sign when z is changed to — z. Therefore, only quadratic con- 
structions of the form (a-b) 2 and (j • b) 2 could, at least in 
principle, enter into the feedback of ijt- 
In an isotropic periodic box we have 

(j-b) = fc f 2 (a-b), (16) 

where fcf can be defined by this relation as the typical 
wavenumber of the fluctuating field. Secondly, we use the re- 
lation (Pouquet, Frisch & Leorat 1976) 

a = -§r(w • u) + ir(j • b)/p (17) 
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for the residual a-effect. This relation describes a fundamen- 
tal form of a-quenching, but there could still be additional 
feedback onto (u) ■ u) and (u 2 ) (Rogachevskii & Kleeorin 
2001, Kleeorin et al. 2002), which is ignored here. With 
these relations, the equation for a becomes explicitly time- 
dependent, 

da 2 /a(B 2 } - ?7t^o(J ■ B) a-a X 

df = ~ 2riwkt 



B 2 



(18) 



where ao = — \t{u ■ u) is the kinematic value of a. Here we 
have expressed the correlation time r in terms of B cq using 
T] t — ^r(u 2 ) and B 2 q = /i po(u 2 ). The full set of equations 
to be solved comprises thus Eqs ( |lT| ) and ([T^). 

A detailed analysis of this set of equations was given by 
Field & Blackman (2002) for the case of the a 2 dynamo and 
by Blackman & Brandenburg (2002) for the aVl dynamo. The 
main conclusion is that for large magnetic Reynolds number 
the large scale magnetic field grows first exponentially such 
that Eq. (^|) is obeyed at all times. This behavior could not 
have been reproduced with an a-effect that is not explicitly 
time-dependent, such as for example Eq. ([[]). The exponential 
growth is then followed by a resistively slow saturation phase, 
just like in the simulations of B01. 

The reason there is this slow saturation phase is that 
Eq. (||) is incompatible with a steady state, where the right 
hand side of Eq. ( |To| ) must vanish, i.e. 

(J • B) = (J • B) + (j • b) = (steady state). (19) 

Since the large scale field is helical, Eq. o) applies and we 
have fci(B ) = fcf(b 2 }, so the large scale field exceeds the 
small scale field by a factor k{ /k\. In the simulations of B01, 
(B")/(b 2 ) was about 4, which is indeed close to kf/k\=5. 
In the beginning of the nonlinear regime, Eq. M) predicts, 

o ^ — 

instead, that (B")/(b 2 ) equals k\/kf, which was 25 times 
smaller in B01. Figure g shows this quantitatively by solving 
Eqs Q and (p|); see Blackman & Brandenburg (2002). 

In order to bring the field ratio from 1/5 to 5 we have to 
remove small scale magnetic helicity resistively. The ques- 
tion is of course what happens if one considers the effects 
of boundaries (both at the equator and at the outer surface): 
can boundaries remove small scale magnetic helicity so that 
the large scale field can saturate at a higher level? This pos- 
sibility was first brought up by Blackman & Field (2000) and 
Kleeorin et al. (2000). 

5. From closed to open boxes 

When the restriction to closed or periodic boxes is relaxed, 
there can be a flux of magnetic helicity through the surface, 
so Eq. ( |io| ) has then an additional term, 

JIT 

— = -2 Wo C - Q, (20) 
at 

where H and C are magnetic and current helicities, respec- 
tively, and Q is the surface integrated magnetic helicity flux. 
In the presence of open boundaries, however, H and Q are no 
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Fig. 2. Evolution of (B } and (b 2 ) (solid and dashed lines, 
respectively) in a doubly-logarithmic plot for an a 2 dynamo 
with ?/t = const for a case with fcf/fci = 10. Note the abrupt 
initial saturation after the end of the kinematic exponential 

2 

growth phase with (B } / (b 2 ) ~ 0.1, followed by a slow sat- 
uration phase during which the field increases to its final su- 
perequipartition value with (B }/(b 2 ) ~ 10. (Adapted from 
Blackman & Brandenburg 2002.) 

longer invariant under the gauge transform A — > A' + V0. 
We use therefore the relative magnetic helicity (Berger & 
Field 1984, Finn & Antonsen 1985), 

H= I (A + Ap)-(B-B P )d^, (21) 
Jv 

where Bp is a potential field used as reference field within 
the volume V, and Ap is its vector potential. Both A and 
Ap can have different (arbitrary) gauges. 

Following Berger & Field (1984), the reference field 
obeys the boundary condition Bp ■ h = B • n, i.e. the nor- 
mal components of both fields agree on the boundary. In slab 
geometry, however, the horizontally averaged mean field has 
to be treated separately and the corresponding reference field 
is (Brandenburg & Dobler 2001) 

Bp = (B) = const, (22) 
so Ap is just a linear function of z, Ap = —z x Bp. Bran- 
denburg & Dobler (2001) found from their simulations that 
most of the magnetic helicity is lost on large scales, where 
the sign agrees with that of the large scale magnetic helicity. 
This was a bit disappointing, because one would have hoped 
that the loss term in Eq. ( pc| ) might supersede resistive losses 
at small scales. That small scale losses can at least in principle 
enhance the large scale field was shown in subsequent simu- 
lations (Brandenburg, Dobler & Subramanian 2002, hereafter 
BDS; see also Fig. |[). The hope is now that this behavior can 
eventually be demonstrated using more realistic geometries. 

6. From boxes to spheres 

Spherical geometry is necessary to assess more realistically 
the helicity transfer through equator and outer surfaces, and 
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Fig. 3. The magnetic energy spectrum for the closed box 
simulation (solid line) compared with the case where small 
scale magnetic energy is removed every 100 time steps, cor- 
responding to St — 1.3 x (fcfu rm s) -1 . Note that the large 
scale magnetic energy (at k = 1) is enhanced relative to the 
reference run, whilst the small scale energy is (as expected) 
reduced. 

the relative contributions from rotation, shear, a-effect, and 
turbulent magnetic diffusion. In a recent paper, Berger & 
Ruzmaikin (2000) estimated that the overall helicity flux at 
the surface would be around 10 47 Mx 2 per 1 1 year cycle. 
This value was also confirmed by BDS, who computed nu- 
merically solutions of the mean-field dynamo equations in 
spherical geometry. The relative magnetic helicity for an ax- 
isymmetric mean field, B = bcj) + V x (a</>), takes the very 
simple form (BDS) 

H n =2 / abdV (23) 

where N denotes the volume of the northern hemisphere. The 
corresponding integrated helicity transfer through the outer 
surface or the equator is 

Q N = -2(f [(5 + Ux B) x (a4>)] • dS, (24) 

J 8V 

where U = r sin 9 VL<j> is the mean toroidal flow. Note that 
Qn = Qns + <?Eq, where Qns is the contribution from 
the outer surface integral and QEq is that from the equato- 
rial plane. 

It is remarkable that even in the presence of just uniform 
rotation alone there is a magnetic helicity flux. For a decay- 
ing dipolar magnetic field, the magnetic helicity flux through 
the equator (from north to south), or, what is the same, into 
the outer surface of the northern hemispheres, or out of the 
surface at the southern hemisphere, is $ 2 per rotation, where 
$ is the magnetic flux through one hemisphere. 

However, once the field is sustained by a dynamo effect 
at a constant amplitude, the helicity flux must be balanced by 



the electromotive force (averaged over one cycle), so 

Q = 2 (S_B)V - 2 Wo (Jj B) V 
= 2a(B 2 )V -277TMo(J-B)y, 

where ryx = r\ + rj t is the total (microscopic plus turbulent) 
magnetic diffusivity. Assuming that the dynamo is saturated 
by a reduction of the residual helicity, see Eq. (|l7|), ao and 
Qn must have opposite signs. This is because saturation re- 
quires that sgn(ao) = — sgn((j • b)), but steady state of 
Eq. ( |l4| ) requires that — sgn((j • b)) = sgn({£ • B)), and 
sgn ((£ ■ B}) = sgn (Q). This can also be seen from a time 
series of magnetic helicity, and the different terms on the right 
hand side the dH^/dt equation; see Fig. Q 

Thus, if Qn (which denotes only the contributions from 
the large scale field), is to be identified with the observed neg- 
ative magnetic helicity flux found on the solar surface (Berger 
& Ruzmaikin 2000, DeVore 2000, Chae 2000), then we must 
conclude that a is negative on the northern hemisphere. This 
scenario, where the main magnetic helicity flux results from 
the large scales, is consistent with the simulations of BD01. 
On the other hand, if small scale magnetic helicity is lost pref- 
erentially at small scale, and the sign of the small scale mag- 
netic helicity is opposite, then a would in that scenario be 
positive on the northern hemisphere. This would be consis- 
tent with the observed negative sign of current helicity on the 
northern hemisphere (Seehafer 1990, Pevtsov et al. 1995, Bao 
et al. 1999, Pevtsov & Latushko 2000), which is plausibly a 
proxy of small scale magnetic helicity. 

Within the framework of mean-field ail dynamo theory, 
a negative a in the nothern hemisphere would explain the ob- 
served migration of the sunspot belts, so one would not need 
to resort to meridional circulation driving the dynamo wave. 
However, there is as yet no well established mechanism to ex- 
plain a negative a (except perhaps magnetic buoyancy with 
shear; cf. Brandenburg 1998). 

Observations do not seem to be able to tell us which of the 
two scenarios is right, because it is difficult to tell whether the 
observed magnetic helicity flux is from large or small scale 
fields. If the observed magnetic helicity flux is from small 
scale, one might wonder why one cannot see the magnetic 
helicity flux from the large scales. On the other hand, if the 
observed magnetic helicity flux is actually already due to the 
large scales, one might expect to see small scale magnetic 
helicity fluxes at higher resolution in the future. 

7. Conclusions 

Magnetic helicity seems to play a much more prominent role 
than what has been anticipated until recently. It has become 
clear that a must satisfy an explicitly time-dependent equa- 
tion. The dynamical a-quenching theory has significant pre- 
dictive power: it describes the different quenching behaviors 
for helical and nonhelical fields, the value of the magnetic 
Reynolds number is explicitly incorporated, and the magnetic 
helicity equation is satisfied exactly at all times. So far, no 
departures between this theory and the simulations have been 
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1.510 1.515 



1.535 




dH /dt 



1.510 1.515 1.520 1.525 1.530 1.535 
t Vl /R 2 

Fig. 4. Time series of the nondimensional ratio H/2/io (M)R 
compared with magnetic energy (in nondimensional units 
and scaled by 1/20). Magnetic helicity is mostly positive in 
the northern hemisphere. The helicity production, S, from 
a-effect and turbulent magnetic diffusion is mostly positive 
and balanced here by a mostly positive magnetic helicity flux 
(dashed and dotted lines refer to the contributions from the 
angular velocity and the a-effect, respectively). 



found. A major restriction of the theory in its present form is 
however the inability to handle cases with spatially nonuni- 
form a-effect. 
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